Superstructures at low spin - high spin transitions 
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Abstract 

In many transition metal compounds, in particular those containing Fe^+ and Co'^^, there oc- 
cur spin-state transitions between low-spin and high-spin (or intermediate-spin) states. We show 
that typical interactions between similar spin-state ions are short-range repulsion, and long-range 
interaction which can have different sign depending on the elastic anisotropy of the lattice and 
on the direction between respective ions. Due to such character of effective interactions at the 
spin-state transitions there may occur different superstructures — ordered arrangement of differ- 
ent spin-states, which in particular may have the form of stripes. The properties of the system 
TlSr2Co05 for which such a superstructure was recently observed experimentally, are discussed 
from this point of view. 

PACS numbers: 71.70.Ch, 61.72.Yx, 64.60. Cn 
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There exist many transition metal compounds in which, due to a competition between 
crystal field splitting and intra-atomic Hund's rule exchange, different spin states lie close 
in energy. Such is often the situation for Fe^"*" and Co'^+. For these ions the low-spin (LS) 
state (occupation of d-levels tg^e^, S=0), intermediate-spin (IS) state iS'lg^l, ^~^) high- 
spin (HS) state {t^g^'^j S=2) can be stabilised at different conditions. In some cases there 
occur spin-state transitions caused e.g. by change of pressure, temperature or composition. 
There are many examples of such transitions among Fe^^-compounds , in particular 
with organic ligands. Among systems with Co^"*", the best known (and still controversial) is 
the case of LaCoOq. For a long time it was assumed that there occurs in it the LS - HS 
,..a,.,tio„ witK ,„c.eas,„g te,„pe..at.e g tat ,n„.e recent ealculat.on. fl and e.pen.e„ts 
^ point rather to the LS - IS transition. 

Usually spin-state transitions occur between homogeneous states having predominantly 
one particular spin state; often they are of first order. However as we shall show below, 
a quite different and interesting situation may occur: the interaction between ions with 
different spin state may be such that at the spin-state transitions an intermediate state 
may be stabilised, in which there appears a spin-state superstructure - an ordered arrange- 
ment of different spin states of the same ionic species, e.g. an alternation of HS and LS 
states. Recently the first example of such superstructure was observed in a layered cobaltite 
TlSr2Co05 and its presence was suggested for NdBaCo205.5 0. 

In the present paper we consider the form of interactions determining the character of 
spin-state transitions, suggest a general mechanism which can lead to spin-state superstruc- 
tures, and qualitatively explain the superstructure observed in TlSr2Co05. 

Conventional description of spin-state transitions 0, 0] relies on interactions induced by 
lattice distortions. Ions with different spin states have different ionic radii, e.g. Co^"*" in the 
LS state in Vl-fold coordination has the ionic radius 0.545 A, and in the HS state — 0.61 A. 
Consequently, if we change the spin state of an ion, e.g. by transforming one of the LS Co's 
into a HS state, it amounts to introducing an "impurity" with larger size into a matrix (the 
so called "sphere in the hole" model sl)- This introduces strain in a crystal, and coupling to 
this strain provides the mechanism of interaction between impurities. In a model of isotropic 
elastic continuum generally used in treating LS - HS transitions 0,0], the only interaction 
remaining is that through the "image forces" due to a stress-free surface of the sample lol |. 
which is attractive. With such an attractive interaction one indeed obtains homogeneous 
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structures, transitions between which may be of first order. 

In real crystals, however, the elastic interactions are more complicated. Thus e.g. in 
weakly anisotropic cubic crystals there appears an interaction between impurities which 
decays rather slowly, as and, most important, which has a different sign in different 

mm-. 



directions 



T(n) 

V{f,f') = -CdQ,Q2-^ (1) 
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where C is a constant of the order of unity, f — f' = \f — f'\ ■ n, Qi and Q2 are the 
"strengths" of impurities {Qi ~ {vi — vq) where Vi is the volume of the impurity and Vq is 
the corresponding volume of the matrix) , and the angular dependence of the interaction (Q) 
is determined by a function of the direction cosines of the vector R = ri — r2'- 

m=nl + nl + nl-l (2) 

Here the elastic anisotropy is characterised by the parameter d = Cu — Cu — 2C44, where 
Cij are the corresponding elastic moduli. We see that e.g. for d < the interaction is 
positive (repulsive) along the [100], [010] and [001] directions and attractive along face and 
body diagonals [110] and [111], and vice versa for d > 0. Thus we see that, in contrast 
to usual assumptions, the interaction between similar states (e.g. between HS ions in a LS 
matrix) may be repulsive, at least in certain directions. 

One can also show that the nearest neighbour (nn) interaction via short-wavelengths or 
optical phonons is typically also repulsive. Thus e.g. in a perovskite lattice, when one puts 
a large HS ion at a certain site, it pushes apart the surrounding oxygens, so that it is more 
favourable to have small LS ions at the nn sites. Consequently, the effective interaction 
between similar ions at the nn sites is repulsive (or equivalently, it is attractive between LS 
and HS states). As is clear from these arguments, one can expect in this case that not only 
homogeneous states (LS or HS), but also certain states with ordered LS and HS ions can 
appear at certain conditions. 

A convenient way to describe spin-state transitions and eventual superstructures is to 
map this situation to an effective lattice gas, or Ising model. Let us consider the situation 



or, 

in the case of TlSr2Co05, IS and HS states 6]; for simplicity we speak below about LS and 



in which each ion can be in two states, e.g. LS and HS (it may be also LS and IS 
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HS). We introduce pseudospin operators cTj = ±1, so that (Xj = — 1 corresponds to a LS state, 
and (Tj = 1 to a HS state at site i. According to our general discussion we can model the 



situation by the effective Ising-type interaction Jij'^i'^ji containing in general short-range 
repulsion, and a long-range part which decays as and which may have different signs 
in different directions; this longer-range interactions depend on the details of the crystal 
structure, elastic anisotropy etc., see Eq.(^. The relative energies of different spin-states, 
e.g. HS vs LS, will be described in this language by the effective "magnetic field" — the 
term /icXj. 

It is known that the Ising model with long-range interactions can give rise to a variety 
of different ordered structures (see the well-known ANNNI model jl^ or the treatment of 
the 2d Ising model with the "Coulomb" interaction in ^^). In this paper we consider a 
simplified model, which nevertheless contains the essential physics, keeping only a small 
number of pair interactions. In contrast to the treatment of jlli ll^l. where one has studied 
systems with fixed concentration of particles (in our mapping — fixed "magnetisation"), in 
our present problem the relative number of different spins (LS and HS ions) is not fixed. 
Consequently we should consider our system not for fixed density ( "magnetisation" ) , but for 
fixed chemical potential ( 'Wgnetic field' ); the .ole of the te„>pe.atu.e can be al.o mapped 
onto a magnetic field jS[. 

We consider a 2d square lattice, modelling the situation in TlSr2Co05, which is the layered 
compound with perovskite-like Co02-plane similar to Cu02-plane in high-T^ cuprates. The 
Hamiltonian of our model is 



H = Ji aiCj + J2 ^i^j - h^(Ti (3) 

{ij)=nn {ij)=nnn i 

Here the first term describes the nearest neighbour interaction, which, as we argued 
above, is repulsive, Ji > 0; the second term is the next nearest neighbour interaction along 
X and y-directions ([100] and [010]) which we also take to be repulsive. The last term in 
Eq.Q describes the difference of the on-site energies of LS and HS states, which can be 
changed e.g. by pressure etc. (For Co^+ the effective field h will be equal to h = QJh — 2A, 
where Jh is the Hund's rule exchange and A is the crystal-field splitting between t2g and 
Cg-levels) . 

It is straightforward to see that for only nn repulsion (Ji > 0, J2 = 0) three states can 
be realized at T = for different values of h: the state with all cTj = —1 (LS-states) and 
with the energy (per site) Ei^s = + h where z is the number of nearest neighbours (in 
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our case z = 4); the state with alternating spins +1,— 1,+1, . . . , forming a two-sublattice 
"antiferromagnetic" structure in a bipartite lattice (we consider below only such a case), 
-E'ls/hs = cind the homogeneous HS state with all cTj = 1 and E^s = ^ — h. Conse- 

quently, we would have jump-like phase transitions between these states with increasing h, 
from the LS state a.t h < —Jiz to an ordered array LS/HS/LS/HS ... for —Jiz < h < Jiz, 
and to a HS state for h > Jiz. Thus we see that this model quite naturally leads to the 
formation of a state with a superstructure of LS and HS states, similar to the one observed 
in 0. 

One important difference of this simple case and the experimental situation observed in 
0] is that the ratio of LS to HS states (IS to HS in the real case of TlSr2Co05) is not 1:1 as 
above, but 1:2, IS states forming diagonal stripes in the 2d square lattice. But just such a 
state appears when one takes into account the second term in the Hamiltonian Eq.Q. 

In order to determine the ground state phase diagram of the model Q we employed a 
Metropolis Monte Carlo al gori thm combined with a single spin flip dynamics to cool the 
states to zero temperature We further compared the energies of all possible periodic 
states with unit cells up to size 5x5. In this way we found four regions with different 
ordered ground states: a ferromagnetic and an antiferromagnetic phase, a 2 x 2 checkerboard 
structure which is degenerate with (2, 2)-stripes, and a (2, 1) stripe phase. The resulting 
phase diagram for h > is shown in Fig.l (as is clear from the form of the Hamiltonian Q, 
the phase diagram for /i < can be obtained by changing a to —a, i.e. by reflecting the 
phase diagram in Fig.l relative to the x-axis and changing + ^ — ). 

Since all the ground states are simple diagonal stripes (except for the checkerboard phase 
which is however degenerate with the (2, 2) stripe phase) we mapped the model onto a one- 
dimensional effective model, which is simpler to analyse and which served us for a detailed 
analysis (see appendix). This study confirms the conclusions presented above and reproduces 
the phase diagram of Fig.l. 

To summarize, away from the boundaries all systematic checks revealed no phases but 
the ones shown in Fig.l. Note however, that on the boundary lines interesting degeneracies 
occur. Thus e.g. on the line h = 4Ji + AJ2 separating the ferromagnetic and the (2,1) 
phase all diagonal stripes (n,l) consisting of n spins a = + and one spin a = —1 are 
degenerate. We stress however that these degeneracies are of no phenomenological relevance 
for the superstructures studied here, since they are strictly confined to the "one-dimensional" 
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phase boundaries. 

The phase diagram obtained above resembles that of the well-known ANNNI model, cf. 
Fig.l in Ref.jlJ], with the difference that because of the ferromagnetic coupling between nn 
spins in one direction in the ANNNI model, (2, 1) stripes in the latter are vertical and not 
diagonal, as in our case (see the appendix for details). Also, our result agrees with Ref . [l5 1 
where the model Eq.Q is discussed for /i = 0. 

From the results presented above we see that in a large part of the phase diagram we 
obtain the phase with (2,1) diagonal stripes; this phase exactly corresponds to the spin-state 
superstructure observed experimentally in TlSr2Co05 j^. 

We also considered the generalized model, including the interaction between the sites 
along diagonals in the plaquette, which according to Eqs. (jll2|) should be attractive (ferro- 
magnetic): the corresponding model in the spin language is 

H = Ji (^i(^j + J2 cncTj + Jd - hY^^i- (4) 

<ij>=nn <ij>=nnn <ij>=diag i 

where Ji, J2 > 0, Jd < 0. 

We found no new ground state phases emerging for this case, but the effect of the diag- 
onal coupling Jrf is merely to shift the phase boundaries in Fig.l by Jd/2 to the right (see 
appendix). 

For a longer range interaction, e.g. the one of the type (0), certain other states may 
appear, depending on the ratio of the constants of the Hamiltonian. The description of 
these states is a formidable problem even in Id-case, especially at finite temperatures 
(see also the results of numerical calculations in |i3l]). We will not discuss all the details 
of this problem here; our main aim is just to demonstrate the possibility of the appearance 
of states with spin-state superstructures at corresponding transitions. We see that such 
superstructures indeed occur quite naturally if we take into account realistic interactions of 
HS and LS states via lattice deformation, including elastic anisotropy and interaction with 
optical phonons. In particular, the superstructure observed experimentally in TlSr2Co05 is 
stable for a certain range of parameters if we include nn and nnn interaction along [100] and 
[010]-directions. 

Concerning the details of the properties of TlSr2Co05 and of the theoretical description 
thereof, a few extra points should be mentioned. We discussed above only the origin and 
type of ordering in the ground state, assuming localised electrons. In reality TlSr2Co05 



undergoes the first order insulator - metal transition with increasing temperature at about 
300 K. We do not discuss this transition in this paper, but one can argue that the energy 
gap in the low-temperature phase may be connected with the occurrence of the spin state 
superstructure. From this point of view one probably may consider this state as a result of 
the formation of the spin-state density wave (SSDW), starting from the high-temperature 
metallic phase. (The treatment presented above would correspond to a strong-coupling limit 
of this picture.) 

Another point worth mentioning is the suggestion j3l that an orbital degeneracy may 
play a role in the properties of TlSr2Co05. In general, indeed, there may exist e^-degeneracy 
in the IS state and tog-degeneracy in both the IS and HS states of Co^+. However from the 
experimental data p] it follows that this orbital degeneracy is predominantly lifted due to 
a rather strong tetragonal elongation of the CoOg-octahedra which exists already in the 
high-temperature phase and which most probably is connected with the layered structure 
of this compound. Indeed, the IS state for this distortion is nondegenerate, and only in 
the HS state there remains a double degeneracy in the t2g-sublevels. But the latter in Co 
compounds is usually lifted by the spin-orbit coupling, which presumably would be the case 
here as well (Jahn- Teller effect by itself would stabilize for the HS state not an elongation, 
but a compression of CoOe-octahedra). This question, however, deserves further study. 

In conclusion, we considered in this paper the properties of systems undergoing spin-state 
transitions, like the ones often observed in materials containing Fe^^ and Co^'^. We argued 
that the effective interaction governing the behavior of such systems - the interaction via 
lattice distortions - is more complicated than the usually assumed attraction between simi- 
lar spin-state ions: typically the nearest neighbour interaction is repulsive, and more distant 
interactions may be of either sign, depending on the elastic anisotropy of the crystal and on 
the direction between respective ions. As a result of such form of interaction, different su- 
perstructures, consisting of ordered distributions of different spin-states, can naturally occur 
in this case. We discussed from this point of view the spin-state superstructure observed in 
the low-temperature phase of TlSr2Co05 and have shown that this superstructure can 
be explained by our model for a certain range of parameters. It would be interesting to look 
for similar superstructures in other systems with comparable energies of different spin-states 
and with spin-state transitions. 
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APPENDIX 



To gain further insight and to corroborate the picture obtained by Monte Carlo calcula- 
tions and by explicitly comparing the energies of different states, we took advantage of the 
in principal one- dimensional structure of the relevant configurations (i.e. diagonal stripes, 
see Fig.l) and mapped the original Hamiltonian onto a one-dimensional Ising model in 
a magnetic field with the Hamiltonian: 

-f^ldim = '^{2JiSiSi+i + (2J2 + Jd)SiSi+2 - hSi + Jd]- (5) 

i 

Here Sj = ±1 represents the normalized sums of the spins on the diagonals and the one- 
dimensional interaction is orthogonal to the stripes. Note that this Hamiltonian contains 
all possible diagonal stripe configurations. Using this one-dimensional model we checked 
systematically for possible diagonal stripe configurations up to unit cell of 18. 

From Eq|3 one sees that — once diagonal stripes are confirmed as ground states — the 
role of Jd is a mere modification of the next nearest neighbour coupling (plus an irrelevant 
constant) and hence a shift of the phase boundaries by Jd/2. 

EqElis also useful to display the similarities to the two-dimensional ANNNI-model 

-f^ANNNI = —'^^{,JlSx,ySx+l,y + J2Sx,ySx+2,y + JoSx,ySx,y+l + HSx,y) (6) 

x,y 

with Ji, J2 < and Jq > as discussed in Ref.^]. Due to the ferromagnetic interaction Jq 
in y-direction, the stripes are vertically oriented in the case of EqlH] and the corresponding 
one-dimensional Hamiltonian is 

-f^ldim = - ^{JlSiSi+i + J2SiSi+2 + HSi + Jq}. (7) 
i 

Hidim Eq.((Zj) differs from i^idim Eq. ^ for = by a factor 2 in the couplings and by 
a constant. The phase diagrams of the models can thus be mapped onto one another by a 
simple rescaling of the axis. 
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FIG. 1: Phase diagram of the model The type of pseudospin ordering is shown for each 
possible phase. (The sign + corresponds e.g. to a high-spin state, and the sign — to a low-spin or 
intermediate-spin state.) Note, that we fixed the energy scale by putting Ji = 1. 
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